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Abstract. Broadband multimode squeezers constitute a powerful quantum 
resource with promising potential for different applications in quantum 
information technologies such as information coding in quantum communication 
networks or quantum simulations in higher dimensional systems. However, the 
characterization of a large array of squeezers that coexist in a single spatial mode is 
challenging. In this paper we address this problem and propose a straightforward 
method to determine the number of squeezers and their respective squeezing 
strengths by using broadband multimode correlation function measurements. 
These measurements employ the large detection windows of state of the art 
avalanche photodiodcs to simultaneously probe the full Hilbert space of the 
generated state, which enables us to benchmark the squeezed states. Moreover, 
due to the structure of correlation functions, our measurements are not affected 
by losses. This is a significant advantage, since detectors with low efficiencies 
are sufficient. Our approach is less costly than tomographic methods relying 
on multimode homodyne detection which is based on much more demanding 
measurement and analysis tools and appear to be impractical for large Hilbert 
spaces. 
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1. Introduction 

The study of correlation functions has a long history and lies at the heart of coherence 
theory p] . Intensity correlation measurements were first performed by Hanbury Brown 
and Twiss in the context of classical optics [2J. Since then correlation functions have 
become an standard tool in quantum optical experiments to study the properties of 
laser beams j3J, parametric downconversion sources [H [5] or heralded single-photons 
[71 [5] . Current state of the art experiments are able to measure correlation functions 
up to the eighth order [5], giving access to diverse characteristics of photonic states. 
The normalized second-order correlation function g^> (0) probes whether the generated 
photons are bunched or anti-bunched, with </ 2 )(0) < 1 being a genuine sign of non- 
classicality |10j . The measurement of all unnormalized moments of a given 

optical quantum state provide complete access to the photon-number distribution 
for arbitrary single- mode input states pQ. Moreover, it is possible to perform a full 
state-tomography with the help of correlation function measurements |1 1 j . 

The measurement of these correlation functions is, in general, performed in a 
time resolved manner g( n '(ti,t2, ■ ■ - t n ). Limited time resolution has been considered 
as a detrimental effect and treated as experimental imperfection [6j. In contrast to 
previous work, we employ the finite time resolution of photo-detectors to gain access 
to the spectral character of broadband multimode quantum states. Our scheme of 
measuring broadband multimode correlation functions of pulsed quantum light is 
especially useful for probing squeezed states. These states are commonly generated via 
the interaction of light with a crystal exhibiting a x( 2 )-nonlincarity, a process referred 
to as parametric downconversion (PDC)[T2j EH HH [15j [16j Hi] or with optical fibers 
featuring a x^-nonlinearity called four-wave-mixing (FWM) p~8| [19]. 

In general the generated squeezed states exhibit multimode characteristics in the 
spectral degree of freedom, i.e. a set of independent squeezed states is created with 
each squeezer residing in its own Hilbert space. This inherent multimode character 
renders these states powerful for coding quantum information, yet the same feature 
impedes a proper experimental characterization in a straightforward manner. Due to 
the sheer vastness of the corresponding Hilbert space, standard quantum tomography 
methods become time-consuming and ineffective. It is neither easy to determine the 
degree of squeezing in each mode, nor the amount of generated independent squeezers. 
Nonetheless, these are the key benchmarks defining the potential of a source for 
quantum information and quantum cryptography applications. In the following we 
investigate how to overcome these issues and elaborate on an alternative approach to 
determine the properties of multimode squeezed states based on measuring broadband 
multimode correlation functions. 

This paper is structured as follows: In section [2] we revisit the general structure 
of multimode twin-beam squeezers drawing special attention — but not restricting 
— to states generated by parametric downconversion and four-wave-mixing. Section 
[3] presents the formalism of correlation functions, introduces the intricacies of finite 
time resolution and defines broadband multimode correlation measurements. Section^ 
combines the findings of sections [2] and [3j We analyze the relation between the number 
of generated squeezers, their respective squeezing strengths and broadband multimode 
correlation functions, which leads us to proposing our scheme for characterizing 
multimode squeezing with the aid of broadband multimode correlation functions. 
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2. Multimode Squeezers 

In a squeezed state of light one quadrature of the field exhibits an uncertainty 
below the standard quantum level at the expense of an increased variance in the 
conjugate quadrature, such that the Heisenberg's uncertainty relation holds at its 
minimum attainable value. The standard description of squeezed states usually 
considers two different types of squeezers: single-beam squeezers and twin-beam 
squeezers. Single-beam squeezers create the squeezing into a single optical mode 
S = cxp (— £a^ 2 + £*a 2 ) , whereas twin-beam squeezers consist of two beams with 

inter-beam squeezing S ab = exp ^— (a^P + (*abj [20] . In these equations £ labels the 

squeezing strength and the operators a\ b' create photons in distinct optical modes. 

In this section we go beyond the standard description and discuss the theory of 
squeezed states, which are generated by the interaction of ultrafast pump pulses with 
nonlinear crystals or optical fibers. Here, we concentrate on the spectral structure of 
the broadband output beams. In general the utilized optical processes, typically called 
optical parametric amplification (OPA) or parametric downconversion (PDC) do not 
generate one but a variety of different squeezers in multiple frequency modes. A whole 
set of independent squeezed beams is generated in broadband orthogonal spectral 
modes within an optical beam. We refer to these states as frequency multimode 
single- or twin-beam squeezers [14j . Here the multimode prefix indicates that more 
than one squeezer is present in the optical beam and the term single- or twin-beam 
identifies whether one squeezed beam or two entangled squeezed beams are created. 
Due to the single-pass configuration of our sources losses are negligible, hence we 
restrict ourselves to the analysis of pure squeezed states. 



2.1. Multimode twin-beam squeezers 

The subject of our analysis is twin-beam squeezing generated by the propagation 
of an ultrafast pump pulse through a nonlinear medium (single-beam squeezers are 
discussed in |Appcndix B[ ) . For simplicity we focus on the collinear propagation of all 
involved fields each generated into a single spatial mode. This description is rigorously 
fulfilled for PDC in waveguides [HJ H2] j but can also be applied for other experimental 
configurations, since the approximation carries all the complexities of the multimode 
propagation in the spectral degree of freedom. If the pump field is undcplctcd, we 
can neglect its quantum fluctuations and describe this OPA process by the effective 
quadratic Hamiltonian (see | Appendix A| for a detailed derivation) 

Hqpa=A \ duj s / (kJif(u a ,Wi)al(u s )aUu)i) + h.c., (1) 



in which the constant A denotes the overall efficiency of the OPA, the function 
f(u> s , u>i) describes the normalized output spectrum of the downconverted beam, which 
— in many cases — is close to a two-dimensional Gaussian distribution. The operators 
a\(u} s ) and a\{uji) are the photon creation operators in the different twin-beam arms, 
in general labelled signal and idler, respectively. 

The unitary transformation generated by the effective OPA Hamiltonian in 
equation ([T]) can be written in the form 



Uqpa = cxp 



(a J du> s j duJif{uj s: uJi)a\{oj s )a\{oJi) + h.c\ 



(2) 
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By virtue of the singular-value-decomposition theorem |23j we decompose the two 
terms in the exponential of equation as 

- rAf(u; s ,u)i) = ~Y] r k ipl(u s )4>l(wi), and 
fc 

_ A* f*(u s ,Wi) = - y^jk4>k{u s )4>k{ui)- (3) 
fc 

Here both {j|)t(w s )} and {4> k (uJi)} each form a complete set of orthonormal functions. 
The amplitudes of the generated modes ?p k (uj s ) and (f> k (uji) are given by the r k £ R + 
distribution. Employing equation ([3]) and introducing a new broadband mode basis 
[24"] for the generated state as: 



A k = J duj s ip k (uj s )a s (uj s ) and B k = / du^fc^a^Wj), (4) 
we obtain the unitary transformation |13) 



Uqpa = cxp 



(g)cxp [r k A{Bl 



h.c. 



= Q$Sf(-r k ). (5) 

fc 

In total the OPA generates a tensor product of distinct broadband twin-beam 
squeezers as denned in [20] with squeezing amplitudes r k , related to the available 
amount of squeezing via: squeezing[dB] = — 10 log 10 (e~ 2rk ) . The Heisenberg 
representation of the multimode twin-beam squeezers is given by independent input- 
output relations for each broadband beam 

i fe =*> cosh(r fe )i fc + sinh(r k )Bl 

B k =>• coBh(r k )B k + sinh(r fc )i|,. (6) 

Note that the squeezer distribution r k and basis modes A k and B k are unique 
and well-defined properties of the generated twin-beam. Their exact form is given 
by the Schmidt decomposition of the joint spectral amplitude — j-Af(uj s , Wj). This 
mathematical transformation directly yields the physical shape of the generated optical 
modes ip k {w s ), 4> k {uJi) with each pair A k and B k being strictly correlated. 

In figure Q] we illustrated one possible squeezer distribution and corresponding 
broadband modes. The joint spectral distribution f(u> s ,u>i) of the generated twin- 
beams shown in figure [T] defines the shape of the broadband signal and idler modes 
A k and B k . In the special case of a Gaussian spectral distribution the form of the 
squeezing modes resembles the Hermite functions. The number of different squeezer 
modes is closely connected to the frequency correlations between the signal and 
idler beam. In the presented case the spectrally correlated beams lead to over 20 
independent squeezers. The total amount of squeezing depends on the constant A 
appearing in the Hamiltonian in equation ([1]) , which is directly related to the applied 
pump power I and the strength of the nonlinearity %( 2 ' in the medium (A cx \/7, x'' 2 '')- 

The OPA state is mainly characterized by the number of squeezed modes and 
the overall gain of the process, both being determined by the distribution of the 
individual squeezing amplitudes r k . In order to analyze the number of generated 
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Figure 1. Visualization of the singular value decomposition in equation The 
frequency distribution — i Af(oj s ,^>i) of the generated state defines the shape of 
the signal and idler modes tp^^uis), (fiki^i) an d the squeezer distribution r^. 



squeezers independently from the amount of squeezing, we split the distribution of 
squeezing weights into a normalized distribution Afc (^ fc \\ = l) that characterizes 
the probability for occupation of different squeezers in the respective optical quantum 
state, and an overall gain of the process B £ R + , quantifying the total amount of 
generated squeezing according to 

r k = B\ k . (7) 

The characterization of these two fundamental properties of a multimode twin-beam 
state is a major experimental challenge. While these states are easily generated in 
the lab, a tomography by means of homodyne detection would require to match for 
each squeezed mode and Bk different local oscillator beams with adapted temporal- 
spectral pulse shapes. Multimode homodyning |25j may provide a route to circumvent 
this difficulty, however an experimental implementation still appears challenging. 



3. Correlation functions 

The n-th order (normalized) correlation function g( n \ti,ti, . . . ,t n ) is generally defined 
as a time-dependent function of the electromagnetic field. For quantized electric field 
operators, it can be expressed as QUI HI H3 

(M-\t 1 )...M-\t n )M+)(ti)...M+\t n )) 

<?(") (h, f 2 , . . . , t n ) = A_ / , ( 8 ) 

E(-)(t 1 )E(+)(t 1 )) . . . (E(-){t n )E(+){t n ) 



and it measures the (normalized) n-th order temporal correlations at different points in 
time. Note that this definition of the correlation functions is independent of coupling 
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losses and detection inefficiencies yielding a loss resilient measure [S]. Realistic 
detectors however, suffer from internal jitter and finite gating times. We accommodate 
for these resolution effects by weighting the correlation function with the appropriate 
detection window T(t) of the applied detectors as presented in [5J, and obtain 

*0 . , *£> . to . 
a) | : 1 b) H c) I 1 

— Optical pulse . 

— Detection gate ITLj ^^^^^^^^^^^^H 

LAL_ LA_, tiAi, 

Figure 2. a) perfect time-resolved detection; b) finite detection gate; c) 
broadband detection gate exceeding the pulse duration giving rise to different 
types of correlation measures. 



9 {n \t x ,t 2 ,...,t n ) = 

J dhTih) ...J dt n T(t n ) (&-\h) . . . &-)(t n )E(+\h) . . . E^(t n )\ 
\ L (9) 

/dt 1 T(ti)(^(-)(t 1 )B(+)(t 1 )).../d* n T(t n )(^(-)(t n )^(+)(t n )) 

If the employed photo-detectors exhibit flat detection windows, exceeding the length 
of the investigated pulses (T(t) — > const.), equation (J9j> can be simplified to 

, Jdt 1 ...dt n (m-\t 1 )...E(-\t n )EM(t 1 )...E^(t n )\ 

= ^ L (10) 

/ dti (E(-)(h)E(+) (h))...Jdt n (M~) (t n )E(+) (t n )) 

This theoretical model is adequate for the detection of ultrafast pulses with standard 
avalanche photodetectors. Furthermore, equation (jlOp exhibits the convenient 
property of time independence and represents our generalized broadband multimode 
correlation function. Despite its similarity to the common correlation functions as 
defined in equation ([8|). the broadband multimode correlation function in equation (|10p 
should no longer be considered as a naive general measure of n-th order coherence. 
In figure [2] we illustrate the main difference between the time-integrated and time- 
resolved correlation measurements. 

Equation (fT0| is still not optimal for our studies of squeezed light fields. We 
transform it further by replacing the electric field operators by photon number creation 
and destruction operators (E^- + \t n ) cx a(t n )) and perform a Fourier transform from 
the time domain into the frequency domain (a(t) — J duj a(uj)e~ tLJt ). Equation (jlOp is 
then rewritten as 

(„) J dujj . ..duj n (a^i) . . ■a t (cj„)a(a;i) . ..a{u n )) 
9 J dwi (at(wi)a(wi)) ... j dio n (at( Wn )a(u;„)) 

_ (:(/d^atnaH)":) ^ 
(J dwaf (w)d(w)) 

in which (: • • • :) indicates normal ordering of the enclosed photon creation and 
destruction operators. In addition we adapt the correlation function to the basis 
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of the measured quantum system, i.e. we perform a general basis transform from a(ui) 
to the basis of the measured multimode twin-beam squeezers A k . This results in: 

Equations (fTDj) . (| 1 1 1) and stress the key difference between time-resolved and 
time-integrated correlation function measurements. While time-resolved correlation 
functions probe specific temporal modes, time- integrating detectors directly measure a 
superposition of all the different modes. This specific feature of broadband multimode 
detection is essential for our analysis. The simultaneous measurement of all different 
optical modes gives us direct loss-independent access to the squeezer distribution of 
the probed state. 

3.1. Broadband multimode cross-correlation junctions 

In the previous section we restricted ourselves to intra-beam correlations. To allow for 
measurements of correlations between different beams we extend our analysis. The 
identification of such inter-beam correlations is of special importance in quantum 
optics and quantum information applications, since they quantify the continuous 
variable entanglement between different subsystems, in our case the analyzed optical 
beams. In scction[5]we have already discussed one of the most employed entanglement 
sources: Twin-beam squeezers. These states arc not only entangled in their 
quadratures, but also in their spectral and spatial degrees of freedom |28| . In order 
to probe higher-order cross-correlations between the two different beams |27j . or 
subsystems a and b of order n and m respectively, we generalize equation to 

(n,m) ( Aa) (a) („) . (6) (6) ( 6 V _ 

[M~\i a) ) ■ ■ ■ Ei-\t^)Ei + \t^) . . . $+>(#>) x £(->(4 b >) . . . Ei + \t^)) 

IWiS) • ■ • >< . . . (wiwi) ' (13) 

Taking into account broadband detection windows — exceeding the pulse duration — 
the above formula can be reformulated as 

9 {n ' m) = ^7 77T, \n r 77-, 77T, /m ' -(14) 



(/ df ^ (t) )" (/ dt E[- ] (t)E[ +) (t))' 



Again we perform the same simplifications as in equation (|lll) in section [31 namely 
we replace the electric field operators by photon creation and destruction operators, 
apply the Fourier transform from the time to frequency domain and finally we adapt 
the measurement basis to the given optical state. We find an extended version of 
equations (fTTj) and ([12]) 



(fdutf(w)&(w)) :: f dw tf(w)b(w) 

9 {n ' m) = ± -A : : 6n 1 (15) 



(16) 



(/dwot(a;)a(w)> U dutf(u)b(u) 
E k A{A k Y:: fe fc ^ 
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Further extensions of cross-correlation measurements to systems consisting of more 
than two different beams are possible [T], but are not necessary within the scope of 
this paper. 

4. Probing frequency multimode squeezers via correlation functions 

Using the theoretical description of squeezers as well as the derived broadband 
multimode correlation functions, we now combine the findings of section [2] and [3] 
We establish a connection between the broadband multimode correlation functions 
and the properties of the squeezing, i.e. the mode distribution A& and the optical gain 
B. 

4-1- Probing the number of modes via g^ 2 ' -measurements 

The foremost important property of frequency multimode squeezers is the number 
of independent squeezers in the generated twin-beam state, which is specified by the 
mode distribution A&. In contrast to the optical gain B, which is easily tuned by 
adjusting the pump power the mode distribution Afc is heavily constricted by the 
dispersion in the nonlinear material and hence — in general — not easily adjustable 
[29] . The effective number of modes in multimode twin-beam state is given by the 
Schmidt number or cooperativity parameter K as defined in [301 131) with 

/v=l/^A|.. (17) 

fe 

Under the assumption of an independent uniform squeezer distribution it directly 
reflects the number of occupied modes. The mode number K of a multimode twin- 
Frequency multimode twin-beam squeezer -measurement 

_ o 



Figure 3. Setup to measure of a multimode twin-beam squeezer. 

beam squeezer can be directly accessed by measuring the broadband multimode g^- 
correlation function in the signal or idler arm as depicted in figure [3] This is a result 
of the structure of the second-order correlation function, which — by using (|12|) and 
([6]) — can be expressed as 

g^ = l+ ^ Smh4 ^) 2 , (18) 
E.sinhV,)] 2 

For our further analysis it is useful to distinguish the low gain from the high 
gain regime, corresponding to low and high levels of squeezing. In the low gain 
regime corresponding to biphotonic states typically referred to in the context of PDC 



Signal 
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experiments sinh(rfc) ~ r^ = BX^ and we are able to simplify equation (|18[) to 
,(2) r~ i i Sfc B A ^k) _ 1 , Sfc ^fc _ 1 . \4 



ff w „ i + 



= 1 + ^- (19) 

Consequently the effective number of modes is directly available from the correlation 
function measurement via K — \/{g^ — 1). For a single twin-beam squeezer (K = 1) 
g( 2 ) = 2, whereas for higher numbers of squeezers (K ^> 1) the contributions from the 
term A^; becomes negligible and g^ approaches one. This direct correspondence 
between g^> and the effective number of modes K is presented in figure 0] (a). 

Another way of interpreting equation (|19p is to approach the correlation function 
measurement from the photon-number point of view. The g^-value of a single twin- 
beam squeezer, which exhibits a thermal photon-number distribution, evaluates to 
g( 2 )= 2. If more squeezers are involved the detector cannot distinguish between 
the different thermal distributions, i.e. it measures a convolution of all the different 
thermal photon streams, which gives a Poissonian photon- number distribution |13U32j . 
In fact one can show that the (^^-correlation function in equation (fT5|) is the 
convolution of the second-order moments of each individual squeezer. 

Once more, we stress that the g^ -measurement does not give access to the 
exact distribution of squeezers A^, but to the effective number of modes under the 
assumption that all squeezed states share an identical amount of squeezing. This is 
a rather crude model and does not fit very well to many experimental realizations. 
Fortunately, there is a common class of squeezed states, for which a much more refined 
mode distribution A& is accessible: In the case of a two-dimensional Gaussian joint- 
spectral distribution f(u) s ,u)i), the distribution Afc is thermal Afc = yl — (jfifi , and 
thus it can be characterized by a single distribution parameter \x |33j . The latter can 
be retrieved from a g^ 1 ' -measurement via /i = -\/2/g( 2 ) — 1, as depicted in figure U] 
(b), where we illustrate how the detection of the g^ 1 -function can provide us directly 
with comprehensive knowledge about the underlying spectral mode structure of the 
analyzed state. 

In conclusion we have shown, that by measuring the second-order correlation 
function g^ of a multimode broadband twin-beam state one can probe the 
corresponding distribution of spectral modes A^ . Our method displays the advantage 
that correlation functions can be measured in a very practical way |34) , resulting in an 
approach that is much easier than realizing homodyne measurements, which require 
addressing individual modes. As a side remark we would like to point out that one can 
also determine the effective number of squeezers from the higher moments 9 {n) ,n>2 
similar to the presented approach, yet g^ is already sufficient for our purposes. 

4-2. Probing the optical gain B of a multimode twin-beam squeezer via g^' 1 ' 
measurements 



In section I4~T1 we determined the number of modes in a loss resilient way by measuring 
g( 2 ) for low gains B. Here we investigate the amount of the generated squeezing 
determined by the overall optical gain B. In order to probe this value the setup has 
to be changed to measure the correlation function g 1,1 ' 1 ^ of the generated twin-beam 
squeezer as presented in figure [5] Using equation (|16p and ^ we obtain for g^ 1 ' 1 ^ the 




Figure 4. a) Plot of the effective mode number K as a function of g( 2 ' for various 
effective numbers of modes, b) Visualization of fi as a function of g( 2 ' for different 
thermal squeezer distributions. 
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Figure 5. Schematic setup to measure g^ 1 '*) of a multimode twin-beam squeezer 
generated via PDC. 



J2 k i sinh 2 (r fe ) sinh 2 (n) + J2k sinh 2 (r fe ) cosh 2 (r fe ) 



[E*sinh 2 (r fc )] 

X; fc sinh 4 (r fc ) 



EfcSinh (r fc ) [E fc sinh 2 (r fc ) 



(20) 



l/(n) 



J l "-1 

The relevant characteristics we exploit from this measurement is its dependence on 
both, the number of modes in the system, as given by the g^-function and the mean 
photon number in each arm, which is closely connected to the coupling coefficient B. 
In the low gain regime (sinh(rfc) w r/~), g^ 1 ' 1 ^ simplifies to 

(i,D ~ i I 1 | Sfc K fa) ■ 1 1 / 2 1) 

ff (2)_l > 

Hence, the optical gain is — in the low gain regime — obtained from the g( 1,1 )- 
measurement via the simple relation B rs 1/ yjW. Mode dependencies of the 
coupling value B only occur at high squeezing strengths, where the relation diverges 
from equation (|21[) and takes on a more complicated form. In figure [5] we plot the 
dependence of the overall coupling value B on g^ 1,1 ^ — as presented in equation ((20)) 
— which takes on a high value for small optical gains B but rapidly decreases when 
the high gain regime is approached. 

In total measuring g( 1 ' 1 ) gives direct loss-independent access to the optical gain 
B. This enables a loss tolerant probing of the generated mean photon number which, 
in the low gain regime, is even independent of the underlying mode structure. 

Taking into account the prior knowledge we gained from section 14. li we can now 
ascertain all parameters needed to fully determine the highly complex multimode 
state. The optical gain B defines not only the photon distribution, but quantifies 
the generated twin-beam squeezing, i.e. the available CV-cntanglement in each mode. 
Note that all modes exhibit different entanglement parameters. Depending on the 
state and its respective mode distribution determined by the (/^-measurement all 
the entanglement could be generated in a single spectral mode where it is readily 
available for quantum information experiments or in a multitude of different squeezed 
modes. Note however, that after the state generation process multiple squeezers cannot 
be combined into a single optical mode by using only Gaussian operations, since 
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Figure 6. The optical gain B plotted as a function of g^ 1 ' 1 ' . For small values 
of B the correlation function g^ 1 ' 1 ' takes on a high value, yet rapidly decreases 
when the high gain regime is approached. 



this operation would be equivalent to continuous-variable entanglement distillation 

030511321. 

5. Outlook 

In this paper we focused on the state characterization of ultrafast twin-beam squeezers 
in the time domain and their experimental analysis. The presented approach however 
is not limited to twin-beam squeezers: 

On the one hand, our measurement technique also applies to probe the squeezing 
of ultrafast multimode single-beam squeezers as presented in |Appcndix B| On the 
other hand, our approach is easily adapted to spatial multimode squeezed states 
[551 1551 HO] . These are characterized by measuring correlation functions that are 
broadband in the spatial domain, in a direct analogy to the spectral degree of freedom 
analyzed in this work. 

6. Conclusion 

We elaborated on the generation of multimode squeezed beams and their 
characterization with multimode broadband correlation functions. We expanded the 
formalism of correlation functions by including the effects of finite time resolution. 
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These extended correlation function measurements serve as a versatile tool for the 
characterization of optical quantum states such as twin-beam squeezers. They provide 
a simple, straightforward and loss independent way to investigate the characteristics 
of multimode squeezed states. Our findings are important for the field of efficient 
quantum state characterization and have already proven to be a useful experimental 
tool in the laboratory [4TJ [34] . 
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Appendix A. Multimode twin-beam squeezer generation via nonlinear 
optical processes 

Appendix A.l. Generation of multimode twin-beam squeezers via parametric 
downconversion 

In the process of parametric downconversion squeezed states are generated by the 
interaction of a strong pump field with the ^^ 2 ^-rionlineaxity of a crystal. Regarding 
the generation of twin-beam squeezers the Hamiltonian of the corresponding thrce- 
wave-mixing process is given by [13l [42] [43] : 

Hpdc = J 2 L dz X {2) E p +) (z,t)Ei-\z,t)E ( -\z,t) + h.c. (A.l) 

where we focused on a collinear interaction of all three beams. In equation (|A.1[) L 
labels the length of the medium, the nonlinearity of the crystal, and E p + \z,t), 
Es (z,t), M~\z,t) the pump, the signal and the idler fields. The electric field 
operators used in equation (|A.1|) are defined as follows 

E { -\z,t)=E^(z,t)=C ( du x exp[-t(k x (oj)z + ujt)}ai(oj), (A.2) 



in which we have merged all constants and slowly varying field amplitudes in the 
overall parameter C. In order to simplify the Hamiltonian we treat the strong pump 
field as a classical wave 

E p + \z,t) =$> E p (z,t) = duj p a(u p )exp[i(kp(uj p )z + ui p t)]. (A.3) 



Here a(oj p ) = A p exp [(w p — /i p ) 2 /(2erp)] is the Gaussian pump envelope function 
generated by an ultrafast laser system, consisting of a field amplitude A p , a central 
pump frequency /i p , and a pump width a p . 

The PDC Hamiltonian in equation (|A.1[) generates the following unitary 
transformation: 



U = exp 



/oo 
dt'H PDC (t') 



(A.4) 



In the low downconversion regime we can ignore the time-ordering of the electric field 
operators [T4] [15] and directly evaluate the time integration. This yields a delta- 
function 2'k8{lj s +LUi — lo p ) and hence allows us to perform the integral over the pump 
frequency oj p . Equation IA.4I can be re-expressed as 

U = exp 

xa(u> s + u>i) exp [lAkz] aj(cj s )a|(wj) + h.c.^j , (A.5) 
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in which Afc = k p {ui s + to^) — k s (u s ) — kiiuii) is the so called phase-mismatch and A 1 
accumulates all constants. Finally, we perform the integration over the length of the 
crystal and obtain 



U = exp 



dujia(ujs + Vi)<j)(uj s ,u}i)al(uj s )al((jJi) + h.c. 



(A.6) 



where (j>(u> s ,u>i) = sine (^rO is referred to as the phasematching function. The 
latter combined with the pump distribution a(uj s + Wj) gives the overall frequency 
distribution or joint spectral amplitude f(uj s ,u>i) of the generated state. The final 
unitary squeezing operator of the downconvcrsion process is 



U = exp 



-r I A I duj s 



dujif (w s ,tx>i)a\{u) s )a\(L)i) + h.c. 



H, 



(A.7) 



The sine function appearing in equation I A. 71 can be approximated by a Gaussian 
distribution 

2" 



/ Ak(oj s ,oJi)L 
[u)s, uJi) = sine I | ps exp 



-0.193 



Ak(oj s , u>i)L 



(A- 



With this simplification the joint frequency distribution /(w s ,Wi) takes on the form 
of a two-dimensional Gaussian distribution. Applying this approximation the exact 
squeezer distribution is accessible as presented in section |U 



Appendix A. 2. Generation of multimode twin-beam squeezers via four-wave-mixing 

In a four-wave-mixing (FWM) process two strong pump fields interact with the x' 3 ^- 
nonlincarity of a fiber to create two new electric fields. If the two generated fields are 
distinguishable the Hamiltonian of the process is given by [44| 



H 



FWM — 



dz (z, t)E { p t ] (z, t)E^(z, t)E$~\z, t) + h.c. 



(A.9) 



Again, we assume a collinear interaction of all interacting beams. The electric fields 
for signal, idler and pump are defined in equations (|A.2[) and (|A.3|) . Performing the 
same steps as in | Appendix A.l wc obtain a similar unitary transformation 



V = exp 



-- I A I duj. 



dwj /FWM(w s ,o;i)aJ(w s )d- (a;;) + h.c. 



H, 



(A.10) 



Equation (|A.10[) resembles equation (|A.7[) with the exception of the joint frequency 
distribution /fwm(w s , Wj) which takes on a more complicated shape in comparison to 
the PDC case 

'Ak((J p ,u> s ,u>i)L\ 



/fwm(w s ,Wi) = / dujp a(u} p )a(uj s + uji - uj p ) sine 



By comparing the unitary transformation in equation (|A.6[) and (| A. 10|) it is apparent 
that the two different processes both create the same fundamental quantum state: 
Multimode twin-beam squeezers. 
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Appendix B. Multimode single-beam squeezers 

In the main body of the paper we discussed the characterization of multimode twin- 
beam squeezers. Here we call attention to the fact that the broadband multimode 
correlation function formalism is also applicable to probe multimode single-beam 
squeezed states. 



Appendix B.l. Generation of multimode single-beam squeezers 



Single-beam squeezers are created by PDC and FWM processes similar to the 
twin-beam states. The difference between the twin-beam and single-beam squeezer 
generation is that in the latter the generated beams are emitted into the same optical 
mode, whereas in the former two different optical modes are generated as discussed in 
|Appcndix A| 

The PDC Hamiltonian generating a single-beam squeezer is given by 



H = 



dz x (2) £v; +) 0, t)E { - } 0, t)E^ ] (z,t) + h.c. . 



(B.l) 



Performing the same steps as in the case of twin-beam generation we obtain the unitary 
transformation 



U = exp 



A / duj s / dw,/(w s ,w ! )a t (w s )a t (w ! ) + /i.c. 



(B.2) 



If the joint spectral distribution f(u) s ,u>i) is engineered to be symmetric under 
permutation of signal and idler, the Schmidt decomposition is given by: 

- T A f{u s ,Ui) =^2r k (f>* k (uj s )<t>* k (uj l ) and (B.3) 



: A*f*{u s ,Wi) = - y^jk<t>k{u s )<j>k{ui) 



(B.4) 



Introducing broadband modes we obtain the multimode broadband unitary 
transformation 



U = exp 



J2r k AlA{-h.c. 



L k 



(g)cxp [r k A\A\-h.c. 



>S(-r k ). 



(B.5) 
I Or 



This is exactly the form of a frequency multimode single-beam squeezed state [1 
written in the Hcisenbcrg picture: 

i fc = cosh(r fc )i fe + sinh(r fe )i^ (B.6) 

Single-beam squeezers are — like twin-beam squeezers — widely employed in quantum 
optics experiments [45j |46j . As in the twin-beam squeezer case the same states are 
generated by properly engineered FWM processes. 
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Appendix B.2. Probing frequency multimode single-beam squeezers via correlation 
function measurements 

In order to characterize the generated states we have to determine the optical gain B 
and mode distribution as in the case of multimode twin-beam squeezers (see section 
|4j. Therefore, we adapt the scheme presented in section |4] and probe the correlation 
functions g^ and g^ as sketched in figure [B~T1 For a multimode single-beam squeezer 



b) 

single-beam squeezer ^'-measurement single-beam squeezer g( 3 )-measurement 



Frequency multimode Frequency multimode 

. 1„ I — -~« r-~,. — n\^l_r 




™P r^near agna| j - Pump ^ sigm| { 



M 



Figure Bl. Schematic setup to measure a) and b) of a frequency 

multimode single-beam squeezer. 



they can be written as: 



foi y^i sinh (rh) 1 

9 l) = 1 + 2 ^=^ + ~ T7— and B.7 

[E fc sinh 2 (r fc )] 2 I>inh 2 (r fc ) 



i/<„) 

,6/ 



(3) _ 1 , 6 Efc s inh (»•*) , Efc sinh ( r *0 



[I>inh 2 (r fc )] [£ fc sinh 2 (r fc )]' 
3 , c EfeSinhV,) 



£ fc sinh> fc ) [£ fe sinh 2 (r fc )]' 

In the single-beam case however g( 2 ' does not directly yield the effective number of 
modes K or thermal mode distribution parameter (i as for the multimode twin-beam 
squeezers in equation (|T9|). A joint measurement of g^ and g^ is necessary, as 
sketched in figure IB2I Clearly the effective mode number K and the thermal mode 
distribution \i are given by the slope s of g^ vs. g( 2 K In figure lB3l we plotted the 
explicit dependence of K and /i on the slope s. Surprisingly the functions exhibit 
almost the same shape as in the twin-beam squeezer case. 

In order to obtain the gain of a multimode single-beam squeezer a single g^- 
measurement is sufficient which is sensitive towards the coupling value B as presented 
in figure [B4l (similar to the g^ 1 '^-measurement in the twin-beam squeezer case). In 
the low gain regime it is given via the relation B = l/s/gi 2 ). Again, while describing 
a different system, the shape of the function B[g^] is very similar to the twin-beam 
squeezer case. 

In total the theoretical description and derivation of multimode single-beam 
squeezers is very similar to the mathematics behind multimode twin-beam states. 
These similarities translate to multimode correlations functions which are able to 
probe the generated optical gain B and mode distribution as in the twin-beam 
case. 




Figure B2. <y' 3 ' as a function of </ 2 ) for various multimode single-beam 
squeezers. The effective number of modes and the thermal mode distributions 
parameter fi of a multimode single-beam squeezer are encoded in the slope. 




b) 




slope 



Figure B3. a) Effective mode number K as a function of the slope of g( 3 ' [<?' 2 ']- b) 
Thermal mode distribution (i as a function of the slope of </( 3 ' [</' 2 '] for multimode 
single-beam squeezed states. 
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Figure B4. Optical gain B as a function of for a multimode singlc- 
squcczcd state. 



